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Abstract 


In  the  paper  by  Mattikalli  et  al. [5],  the  stability  of  an  assemblage  of  frictionless  contact¬ 
ing  bodies  with  uniform  gravity  was  considered.  The  problem  of  finding  a  stable  orientation 
for  such  an  assembly  was  formulated  as  a  constrained  maximin  problem.  A  solution  to 
the  maximin  problem  yielded  an  orientation  of  the  assembly  that  was  stable  under  gravity; 
however,  if  no  such  orientation  existed,  then  the  solution  to  the  maximin  problem  yielded  the 
most  stable  orientation  possible  for  the  assembly.  The  maximin  problem  was  solved  using 
a  numerical  iteration  procedure  that  solved  a  linear  program  for  each  step  of  the  iteration. 
In  this  paper,  we  show  that  the  stability  problem  can  be  considered  a  variant  of  standard 
zero-sum  matrix  games.  A  solution  to  the  maximin  problem  can  be  found  by  solving  a 
single  linear  program. 


Figure  1 :  Body  B  is  fixed  in  place,  while  bodies  A  and  C  are  free  to  move,  (a)  The  assembly 
is  unstable  for  the  current  choice  of  g.  (b)  The  assembly  has  been  reoriented  so  that  it  is 
stable. 

1  Introduction 

In  Mattikalli  et  al.[5 ]  we  showed  how  a  stable  orientation  for  an  assembly  could  be  for¬ 
mulated  as  the  solution  to  a  constrained  maximin  problem.  We  define  an  assembly  as  a 
collection  of  frictionless  contacting  rigid  objects,  one  or  more  of  which  is  assumed  to  be 
fixed  in  place  (for  example  a  floor,  a  supporting  surface  such  as  a  table,  or  an  object  held 
by  a  gripper).  All  objects  in  the  assembly  are  initially  motionless  and  are  acted  upon  by 
an  external  force  mg  where  m  is  an  object’s  mass  and  g  is  a  unit  vector  pointing  straight 
down.  If  the  objects  remain  motionless  under  the  influence  of  the  gravity  field,  we  say  the 
assembly  is  stable.  Otherwise,  the  assembly  is  unstable. 

An  assembly  that  is  unstable  in  one  orientation  might  be  made  stable  by  a  change  of 
orientation  (figure  1).  Given  an  unstable  assembly,  we  would  like  to  be  able  to  reorient 
the  assembly  to  make  it  stable,  if  possible.  However,  rather  than  actually  rotate  objects  in 
our  frame  of  reference,  we  will  instead  choose  a  different  unit  gravity  vector  for  which  the 
assembly  is  stable,  if  such  a  gravity  vector  exists  (figure  2).  When  we  say  we  are  searching 
for  a  stable  orientation  then,  we  mean  that  we  are  searching  for  a  gravity  direction  g  for 
which  the  assembly  is  stable.  If  it  turns  out  that  the  assembly  is  unstable  no  matter  what 
direction  of  gravity  we  choose,  then  we  would  like  to  determine  an  orientation  (that  is,  a 
gravity  direction)  which  makes  the  orientation  as  stable  as  possible.  We  will  define  the 
metric  used  to  measure  stability  in  section  5. 

The  formulation  we  use  for  finding  the  most  stable  orientation  is  based  on  energy 
considerations.  If  an  assembly  is  initially  at  rest,  the  kinetic  energy  T  of  the  assembly  is 
zero.  If  the  assembly  is  unstable  and  begins  to  move,  then  T  must  increase.  Since  the  only 


Figure  2:  Rather  than  actually  rotate  objects,  a  new  gravity  direction  g  is  chosen  so  that  the 
assembly  in  figure  la  is  stable. 

external  force  is  gravity,  which  is  conservative,  if  T  increases  the  potential  energy  U  of  the 
assembly  must  decrease.  Thus,  if  every  motion  that  does  not  cause  interpenetration  between 
objects  is  “uphill”  (that  is,  causes  U  to  increase),  then  the  configuration  is  stable.  Our  goal 
then  is  to  find  a  direction  of  gravity  so  that  all  allowable  motions  are  “uphill.”  If  this  is 
not  possible,  we  want  to  orient  gravity  to  minimize  the  steepness  of  the  most  “downhill” 
allowable  motion. 

The  optimal  direction  for  g  in  Mattikalli  etal.[5]  was  defined  in  terms  of  the  solution  to 
a  constrained  maximin  problem.  The  maximin  problem  was  solved  by  numerical  iteration. 
Each  step  of  the  iteration  involved  solving  a  linear  program.  We  will  show  that  a  solution  to 
this  maximin  problem  can  be  found  by  solving  a  single  linear  program,  eliminating  the  need 
for  an  iterative  solution  method.  The  insight  into  this  reduction  lies  in  viewing  the  maximin 
problem  as  a  variant  of  two-person  zero-sum  matrix  games.  Two-person  zero-sum  matrix 
games  were  first  shown  to  be  equivalent  to  linear  programs  by  Dantzig.  Mattikalli  et  al. 
discusses  previous  work  on  stability  problems. 

2  Motion  Constraints 

We  will  represent  possible  motions  of  an  assembly  in  terms  of  virtual  displacements.  Let 
6p,-  =  (6iv,  60,)  represent  a  displacement  of  the  ith  body  in  the  system,  with  6r,  and  60, 
vectors  in  R3.  The  vector  6r,  denotes  a  translational  displacement  of  the  ith  part,  while  60/ 
denotes  a  rotation  of  magnitude  ||60,j|  of  the  body  around  its  center  of  mass.  The  axis  of 
the  rotation  is  along  the  60,  direction. 

Because  of  contact  between  bodies,  not  all  displacements  are  allowable.  Consider 


Figure  3:  Contact  between  bodies  A  and  B.  Motion  constraints  are  formulated  in  terms  of 
the  relative  motion  of  the  bodies  at  points  d  and  d'. 

figure  3  where  bodies  A  and  B  contact.  If  body  A  undergoes  a  displacement  6pa  =  (6ra,66a), 
then  point  d,  as  attached  to  body  A,  undergoes  a  particular  displacement  <5da.  Similarly,  a 
displacement  6pb  of  body  B  causes  a  displacement  <5d*  of  point  d,  as  attached  to  body  B.  To 
prevent  interpenetration  from  occurring,  the  relative  displacement  <5da-<5dfc  cannot  have  any 
component  opposite  the  unit  normal  direction  n.  We  can  express  this  as  the  constraint 

n  •  (6da  -  6db)  >  0.  (1) 

Similarly,  we  also  need  to  prevent  interpenetration  from  occurring  at  point  d'  by  requiring 
that  n  •  (6d^  -  MJ,)  >  0. 

Note  that  we  do  not  have  to  generate  a  constraint  for  every  point  of  contact  between 
two  bodies.  It  is  only  necessary  to  place  constraints  on  the  vertices  of  the  convex  hull  of  the 
contact  region.  In  this  paper,  we  assume  that  all  contact  regions  are  polygons  (or  degenerate 
polygons),  which  means  that  we  have  to  formulate  constraints  for  only  finitely  many  contact 
points.  Baraff[2]  discusses  the  difficulties  in  dealing  with  nonpolygonal  contact  regions. 
We  will  assume  that  the  motion  constraints  of  the  assembly  can  be  expressed  by  a  finite 
number  of  constraint  inequalities  in  the  form  of  equation  (1),  all  of  which  must  be  satisfied. 
(Palmer[6]  discusses  the  complexities  that  arise  when  only  a  subset  of  the  motion  constraints 
need  to  be  satisfied  by  a  given  motion.) 

For  the  remainder  of  the  paper  then,  we  will  consider  assemblies  whose  motion  con¬ 
straints  are  expressed  in  terms  of  m  contact  points  between  the  bodies.  Let  the  t'th  contact 
point  of  the  assembly  be  a  contact  between  bodies  A  and  B  at  the  point  d,  in  a  global  frame 
of  reference.  Let  n,  denote  the  unit  surface  normal,  pointing  outwards  from  B  towards  A  at 
d(.  Let  ca  and  c *  denote  the  position  of  the  center  of  mass  of  bodies  A  and  B  respectively. 
If  A  undergoes  a  displacement  6pa  =  (6ra,  S0a )  then  d„  as  attached  to  A,  undergoes  the 
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displacement 


6ra  +  66 a  x  (d,-  -  ca). 

Similarly,  for  a  displacement  6pb  =  (6rb,  60b)  of  body  B,  d,’s  displacement,  as  attached  to 
B,  is 

6rb  4-  60 b  x  (d,-  -  cb). 

The  motion  constraint  at  the  ith  contact  point  is  therefore 

n,  •  ( 6ra  +  60a  x  (d(-  -  ca)  -  6rb  -  66b  x  (d,-  -  c*))  >  0.  (2) 

Note  that  if  contact  occurs  between  a  body  A  and  a  fixed  object  B,  we  write  a  constraint 

without  referring  to  a  displacement  of  the  fixed  body,  as 

n,  •  (6ra  +  60a  x  (d/  -  cfl))  >  0.  (3) 

To  simplify  bookkeeping,  we  do  not  count  fixed  objects  as  bodies  in  our  assembly;  rather, 

we  simply  note  when  regular  movable  objects  are  in  contact  with  fixed  objects,  and  generate 
the  appropriate  motion  constraint,  such  as  equation  (3). 

Since  each  constraint  is  a  linear  inequality  on  the  6r  and  66  variables,  we  can  express 
the  simultaneous  satisfaction  of  all  the  constraints  as  one  large  linear  system.  If  the  vector 
<5p  denotes  the  virtual  displacements  of  an  assembly  with  n  bodies,  that  is 

^  (Sr*  N 
60  { 

6p=  i 

6r„ 

\  60n  t 

then  we  express  all  m  motion  constraints  by  writing 

J<5p  >  0  (4) 

where  J  is  an  m  x  6n  matrix  (since  each  displacement  pair  (<5ry,  66 j)  has  six  components)  and 
0  is  an  appropriately  sized  vector  of  zeroes.  (We  will  denote  row  vectors,  column  vectors, 
and  matrices  whose  entries  are  all  zero  simply  by  0  throughout  this  paper.  The  dimension 
of  0  should  be  clear  from  the  context  in  which  it  occurs.)  The  coefficients  of  J  are  computed 
according  to  the  constraint  equations  (2)  and  (3). 

Using  this  notation,  we  can  say  that  a  legal  motion  for  an  assembly  is  a  displacement 
<5p  that  satisfies  J(5p  >  0.  Note  that  the  displacement  6p  =  0  always  yields  a  legal  motion 
(the  null-motion).  Also,  if  <5p  ^  0  is  a  legal  motion  and  a  is  a  nonnegative  scalar,  then  a<5p 
is  also  a  legal  motion,  but  —  a6p  is  not. 

3  Determining  Stability 

Determining  if  an  assembly  is  stable  under  a  particular  direction  of  gravity  is  fairly  straight¬ 
forward.  In  this  section,  we  describe  two  different  methods  for  determining  if  an  assembly 
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is  stable.  The  first  method  is  based  on  potential  energy  considerations,  while  the  second 
method  considers  the  contact  forces  that  arise  at  contact  points.  Both  methods  involve  linear 
programming.  The  latter  method  was  used  in  work  by  Blum,  Griffith,  and  Neumann[3], 
and  is  not  limited  to  frictionless  assemblies.  In  the  next  section,  we  show  how  both  of 
these  methods  for  determining  stability  can  be  modified  to  find  a  stable  orientation  for  the 
assembly  (if  it  exists). 

3.1  Potential  Energy 

Suppose  that  an  assembly  undergoes  a  virtual  displacement  <5p.  The  change  in  potential 
energy  8U  corresponding  to  the  displacements  <5p,  =  (<Sr/,  69 ,)  of  the  n  bodies  is 


r,  (5) 

f=i 


where  M,  is  the  mass  of 

the  ith  body.  If  we  define  the  matrix  M  as  the  3 

X 

6 n  matrix 

(Mi 

0 

0 

0  0  0.. 

■  Mn 

0 

0 

0 

0 

M  = 

\  o 

Mj 

0 

0  0  0.. 

0 

0 

0 

0 

0  ’ 

l  0 

0 

M, 

0  0  0.. 

0 

0 

Mn 

0 

0 

oj 

then  we  can  write 

n 

EM,<5r,= 

1=1 

M<5p 

(6) 

and  thus 

8U 

=  -g  •  M6p  = 

=  —  grM<5p. 

(7) 

As  stated  in  »ne  introduction,  if  for  a  given  gravity  direction  g  all  legal  motions  yield 
6U  >  0,  then  the  assembly  is  stable.  We  are  therefore  interested  in  knowing  the  minimum 
value  that  8U  can  assume  over  all  legal  motions.  If  we  let  z.  denote  this  minimum,  by 
writing 

z  =  min  -grM<5p  (8) 

Jop>0 

then  the  structure  is  stable  if  z  >  0.  Note  however  that  the  null-motion  Ap  =  0  is  always 
legal,  and  yields  6U  =  0.  Thus,  z  is  bounded  above  by  zero,  and  we  can  say  simply  that  the 
structure  is  stable  if  z  =  0. 

The  value  z  can  be  determined  by  linear  programming.  However,  as  it  stands,  if  z  #  0, 
then  there  must  exist  a  legal  <5p  for  which  6U  <  0.  In  this  case,  the  minimum  value  of  6U 
is  —  oo,  since  a  displacement  ot  a6p  is  legal  and  yields  an  energy  change  of  n6U  <  0  for 
any  a  >  0.  This  is  a  consequence  of  the  constraint  J0p  >  0,  which  constrains  the  motion 
direction  <5p,  but  not  its  magnitude. 

Anticipating  future  development,  it  is  useful  to  bound  the  magnitude  of  the  displace¬ 
ments  <5p  considered  in  equation  (8).  Since  we  would  like  to  be  able  to  use  linear  program¬ 
ming  techniques,  we  would  like  to  bound  6p’s  magnitude  with  linear  constraints.  We  can 
do  this  straightforwardly  by  redefining  z  as 


The  infinity  norm  ||  vjloo  of  a  vector  v  is  the  maximum  absolute  value  overall  the  components 
of  v.  The  condition  ||6p||oo  <  1  constrains  all  components  of  6p  to  have  magnitude  less 
than  one.  Equation  (9)  can  then  be  solved  by  linear  programming.  If  the  solution  z  is  zero, 
then  the  assembly  is  stable.  Otherwise,  z  is  a  (finite)  negative  value,  and  the  assembly 
is  unstable.  (Note  however  that  the  displacement  6p  which  yields  the  minimal  z  =  6U 
only  approximately  indicates  the  direction  of  impending  motion  of  the  assembly.  In  order 
to  exactly  determine  the  impending  motion  direction,  it  is  necessary  to  solve  a  quadratic 
programming  problem!  1  ].) 

3.2  Contact  Forces 

Instead  of  looking  at  motion  directions  which  decrease  potential  energy,  we  can  consider 
the  contact  force  that  arises  at  each  of  the  m  contact  points  of  the  assembly.  Since  we  are 
dealing  with  frictionless  contacts,  we  know  that  the  contact  forces  will  act  normal  to  the 
contact  surfaces.  Thus,  at  the  /th  contact  point,  we  consider  a  contact  force  /,n,  that  acts  on 
body  A  of  the  contact,  and  a  contact  force  —  f, n,  that  acts  on  body  B  of  the  contact,  with  /, 
the  unknown  scalar  magnitude  of  the  force.  Since  is  directed  from  B  towards  A,  and  since 
contact  forces  must  be  repulsive,  the  magnitude  /,  must  be  nonnegative;  that  is,  /,  >  0. 

Let  the  vector  of  contact  force  magnitudes  /;  be  denoted  by  f.  The  net  force  Fj  G 
acting  on  the  jth  body  of  the  assembly  can  be  written  as 

=  +  (,0) 

t=i 

where  Sj,  is  either  1 ,  -  1 ,  or  zero.  If  the  jth  body  is  not  involved  in  the  /th  contact,  then  s„  is 
zero.  If  the  contact  force  exerted  on  the  jth  body  from  the  /th  contact  point  is  f,h,,  then  s,, 
is  1.  Otherwise,  the  contact  force  acting  on  the  jth  body  is  —  /,n„  and  Sj,  is  —  1 . 

The  net  torque  r,  G  R3  acting  on  the  jth  body  of  the  assemb'y  is  similarly  written  as 

n 

Ti  =  E^.(d<  ~c/)  x  M  00 

(=1 

where  d,  is  the  location  of  the  /th  contact  point,  and  c ,  is  the  location  of  the  center  of  mass  of 
the y'th  body.  The  scalars  Sj,  are  the  same  as  in  the  previous  equation.  The  r7  are  independent 
of  g  since  a  uniform  gravity  field  does  not  exert  a  torque.  If  we  define  the  6/i-vectors  Q  and 
G  as  the  collections 


f  F,  ^ 

(  M\ g  ) 

Q  = 

ti 

and  G  = 

0 

F„ 

l  T"  ) 

Mng 

\  °  ; 

we  can  write 

Q  =  A/  +  G 


(12) 


(13) 
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where  A  is  a  6n  x  m  matrix  whose  coefficients  are  given  by  equations  ( 1 0)  and  (11). 

Because  the  assembly  is  frictionless,  its  impending  motion  is  completely  determined[4]. 
If  there  exist  repulsive  contact  forces  such  that  the  net  force  and  torque  on  every  body  is 
zero,  then  such  forces  will  arise  at  the  contact  points,  and  the  assembly  is  stable  and  will 
not  move.  Thus,  we  can  determine  stability  by  simply  checking  whether  there  exists  /  such 
that 

Q  =  A/  +  G  =  0  and  /  >  0.  (14) 

The  existence  of  a  suitable  /  can  Lc  determined  by  linear  programming. 


4  Finding  a  Stable  Orientation 

Suppose  that  for  an  assembly  and  a  given  gravity  direction  we  find  that  the  structure  is 
unstable,  using  either  of  the  two  methods  described  in  the  previous  section.  Can  we  find  a 
new  direction  for  g  which  will  make  the  structure  stable?  In  this  section,  we  show  how  the 
contact-force  formulation  to  determine  stability  can  be  trivially  modified  to  find  a  value  for 
g  (if  it  exists)  which  makes  the  assembly  stable. 

In  section  3.2,  we  were  searching  for  a  vector  /  so  that 

Q  =  A/  +  G  =  0  and  />0  (15) 

where  G  depended  upon  the  known  value  of  g.  If  however  we  treat  g  as  an  additional 
unknown,  all  we  need  to  do  is  simply  check  and  see  if  there  exist  values  for  f  and  g  that 
satisfy  Q  =  0,  indicating  that  the  assembly  is  stable.  That  is,  if  we  can  find  vectors  g  and 
/  such  that 

Q  =  A/  +  G  =  0,  f>  0  and  ||g||2  =  1  (16) 

(where  G  is  defined  in  terms  of  g  by  equation  ( 1 2))  then  the  assembly  is  stable  in  orientation 
g.  If  equation  (16)  has  a  solution,  we  can  find  it  by  linear  programming,  although  a  slight 
modification  is  required. 

The  constraint  that  g  be  a  unit  vector  (that  is,  ||g|[2  =  1)  cannot  be  enforced  in  a  linear 
program.  However,  it  is  not  necessary  to  search  among  only  unit  vectors;  it  is  merely 
necessary  to  make  sure  that  we  search  among  all  possible  directions.  We  can  do  this  by 
considering  vectors  g  such  that  ||g|ji  =  1.  (For  a  vector  v,  |jv||i  =  ]T(  |v,|.)  However,  since 
the  set  of  vectors  ||g||i  =  1  is  nonconvex,  we  need  to  split  it  up  into  pieces.  Intuitively, 
||g||i  =  1  forms  a  “unit  diamond”  about  the  origin,  consisting  of  eight  planar  facets.  Let 
the  set  Si  be  defined  as  the  set  of  vectors  g  =  (gx,  gy,  gz)T  satisfying 

gx,gy,gz>  0  and  gx+gy  +  gz=  1-  07) 

Similarly,  define  Sz  to  consist  of  vectors  satisfying 

-gx,gy,gz>  0  and  ~gx  +  gy  +  gz=  1  (18) 

and  so  on  through  all  the  eight  sign  permutations  of  gx,  gy  and  gz.  Using  this  notation,  we  can 
look  for  a  g  that  makes  the  assembly  stable  by  seeing  if  any  of  the  eight  linear  programs 

Q  =  A/  +  G  =  0,  />  0,  and  g€$(l<i<8)  (19) 
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have  a  solution  /  and  g.1 

Note  that  if  an  assembly  can  be  made  stable,  we  find  the  value  of  g  which  makes  it  stable. 
However,  when  no  such  g  exists,  all  we  know  is  that  the  linear  program  is  unsatisfiable. 
Even  if  an  assembly  has  no  stable  orientation,  it  would  still  be  desirable  to  know  what 
orientation  g  comes  the  closest  to  making  the  assembly  stable.  In  the  next  section,  we  will 
show  how  such  a  direction  g  can  be  found  by  linear  programming  using  the  potential  energy 
formulation. 


5  Finding  the  Most  Stable  Orientation 

Modifying  the  contact-force  formulation  to  find  a  stable  orientation  (if  it  exists)  was  straight¬ 
forward.  In  working  with  the  potential  energy  formulation  though,  we  are  not  limited  to 
simply  finding  a  stable  orientation,  or  reporting  that  the  assembly  is  unstable.  Instead,  we 
can  modify  the  potential  energy  formulation  so  that  we  can  find  either  a  stable  orientation, 
or,  for  unstable  assemblies,  the  most  stable  orientation  possible. 

Let  us  recast  2  as  a  function  of  g  by  writing 

2(g)  =  min  — grM(5p.  (20) 

J<5p>0 

ll«Pl|oo<» 

Ifg  is  an  orientation  for  which  the  assembly  is  stable,  then  z(g)  =  0.  Otherwise,  the 
assembly  is  unstable,  and  2(g)  <  0.  We  will  use  the  function  2(g)  as  a  measure  of  the 
instability  of  an  assembly  under  gravity  g.  Since  there  may  be  no  value  g  for  which  z(g)  =  0, 
we  will  search  for  a  value  of  g  that  maximizes  2(g).  If  this  maximum  is  zero,  then  we 
will  have  found  a  stable  orientation.  Otherwise,  we  will  have  found  the  “most”  stable 
orientation,  as  defined  by  the  metric  2(g).  In  the  remainder  of  this  paper,  we  will  restrict 
g  to  lie  in  S i,  by  writing  g  >  0  and  ||g||i  =  £j=i  =  L  In  searching  for  the  most  stable 
assembly,  we  will  have  to  perform  eight  different  searches;  one  for  each  partition  5,.  All 
statements  and  methods  made  hereafter  involving  g  6  St  can  be  applied  to  the  other  seven 
cases  of  g  e  S;. 

5.1  Maximin 

To  find  the  most  stable  orientation,  we  are  trying  to  solve  a  maximin  problem.  That  is,  we 
are  trying  to  solve 

max  2(g)  =  max  min  -grM6p  .  (21) 

g>o  g>o  l  J£p>0 

Ef! 1=1  S«|  =  I  \l|6plloo<l  / 

1  We  could  also  partition  gravity  by  considering  the  unit  cube  of  directions  ||g|loo  =  1  The  natural  division 
here  would  be  to  use  sets  Si  through  Sf,,  with  S|  defined  by  gx  =  1  and  - 1  <  gy,gz  <  1  and  S’  by  gx  =  - 1 
and  -1  <  gy,gz  <  I,  and  similarly  for  S3  through  S6.  This  would  work  just  as  well  as  considering  the  set 
j|g|li  =  I.  The  only  reason  for  using  the  metric  ||g||i  =  I  in  this  section  is  that  we  are  forced  to  use  this 
metric  in  subsequent  sections. 
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Constrained  maximin  problems  are  in  general  hard  to  solve.  However,  problem  (21) 
has  a  form  similar  to  a  maximin  problem  that  is  solvable  by  linear  programming.  Given 
an  m  x  n  matrix  A,  the  minimax  theorem  of  matrix  games,  first  proved  by  Von  Neumann, 
states  that 


max 

f  T  ) 

min  y  Ax 

=  min 

f  T  ) 

max  y  Ax 

x>0 

,  y>°  / 

y>0 

,  *>0  1 

S  jt,=1 

i  J 

£yi-i 

\Zx,=\  / 

where  x  and  y  are  vectors  in  Rm  and  R”  respectively.  Furthermore,  the  value  of  x  for  which 
the  maximum  on  the  left  is  attained  can  be  found  by  solving  a  linear  program;  the  solution  of 
the  dual  linear  program  gives  the  value  of  y  for  which  the  minimum  on  the  right  is  attained. 

The  maximin  problem  problem  (21)  is  similar  to  problem  (22)  in  that  the  constraints  on 
g  are  g  >  0  and  ]Tg,  =  1 .  However,  the  constraints  on  the  inner  variables  <$p  have  quite  a 
different  form.  It  turns  out  that  a  variant  of  the  linear  program  used  to  solve  problem  (22)  for 
x  can  be  used  to  solve  problem  (21)  for  g.  In  section  6  v/e  will  exhibit  a  pair  of  dual  linear 
programs  which  find  a  solution  g  to  problem  (21)  and  a  solution  6p  to  the  dual  problem 


min 

Jfp>0 

n«wioo<i 


max  -grM«p 
W=i 


It  will  also  turn  out  that 


max  z(g)  =  max 
8>o  g>o 

E«,=  l  £«;=! 


(min  -grM6p 

J«p>0 

Halloo  <  i 


min  max  —  g^Mtfp 
J*p>o  1  g>o  I 

IIWloo  <  I  \Zt!,  =  I  / 


(23) 


(24) 


For  now  however,  it  will  be  more  instructive  to  simply  hope  that  equation  (24)  holds, 
and  use  physical  intuition  to  formulate  a  linear  program  that  solves  equation  (23)  for  <5p. 
Assuming  that  equation  (24)  holds,  the  solution  of  the  dual  of  this  linear  program  will  yield  a 
vector  g  which  maximizes  problem  (2 1 ).  The  intuition  which  allows  us  to  directly  formulate 
the  necessary  linear  program  lies  in  viewing  problem  (23)  as  a  competition,  or  game  (just 
as  problem  (22)  is  viewed  as  what  is  known  as  a  “two-person  zero-sum  game”). 


5.2  A  Particle  versus  Gravity 

For  illustrative  purposes,  let  us  greatly  simplify  the  problem.  Our  assembly,  for  the  moment, 
consists  of  a  single  particle  in  R2,  with  degrees  of  freedom  brx  and  bry  and  unit  mass.  Our 
gravity  vector  likewise  has  two  components  gx  and  gy.  To  begin  with,  we  will  assume  that 
there  are  no  constraints  on  6rx  or  6ry,  except  for  the  bounds  l^l  <  1  and  j<5ry|  <  1.  We 
will  search  for  the  most  stable  orientation  by  finding  the  solution  of 


min 

ISrtUtry\<l 


(max  -L 

K',Ky>0 
K,+Kv=  1 
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where 


L=(t)‘(^)=g>I  +  gAr 

Let  us  simplify  this  by  removing  the  minus  sign;  this  swaps  the  “min”  and  “max”  functions, 
yielding 

(  \ 

max  min  L  (25) 

l*r.l.l«'vl<l  l  s,,£y>0  v  ' 

\Sx+«y=l  / 

Normally,  we  think  of  forces  such  as  gravity  as  having  an  “upsetting”  effect  on  objects 
in  that  they  cause  objects  to  move  unless  balanced  by  other  forces.  Objects,  on  the  other 
hand,  have  inertia — they  seek  to  stay  where  they  are  unless  compelled  to  move  by  forces. 
We  can  view  equation  (25)  as  a  game  by  changing  our  viewpoint  slightly.  We  will  regard 
both  the  force  of  gravity,  and  the  particle,  as  players  in  a  game,  each  with  a  particular 
purpose.  On  one  side,  we  have  player  “Particle.”  Rather  than  stand  still,  Player  Particle 
would  like  to  move  as  much  as  possible  in  the  direction  of  gravity  chosen  by  Player  Gravity. 
Mathematically,  player  Particle  seeks  to  maximize  L.  As  a  result.  Particle  is  very  lazy. 
While  Particle  is  quite  happy  to  make  L  positive  by  moving  downhill  (that  is,  in  the  direction 
chosen  by  Gravity),  Particle  is  so  lazy  that  rather  than  move  uphill  the  slightest  bit,  making 
L  negative.  Particle  would  rather  just  stand  still. 

Conversely,  Player  Gravity’s  goal  is  to  pick  a  gravity  direction  so  that  L  is  as  small  as 
possible.  In  terms  of  the  game,  Gravity  needs  to  pick  a  nonzero  direction  that  causes  the 
particle  to  want  to  move  as  little  as  possible  (which,  as  we  said,  is  opposite  our  normal 
depiction  of  gravity  as  an  upsetting  force  on  systems).  Gravity  is  therefore  annoyed  when 
the  particle  actually  moves  in  the  direction  of  gravity,  since  this  makes  L  negative.  Gravity 
would  prefer  the  particle  not  to  move  at  all,  which  makes  L  zero. 

Now  suppose  that  the  game  is  played  without  any  constraints  on  Particle’s  motion. 
(Clearly,  Particle  has  the  advantage  here.)  Particle  reasons  that  by  not  moving  at  all  (that 
is,  by  choosing  brx  =  bry  =  0),  a  value  of  L  =  0  is  attained.  However,  since  Particle  knows 
that  Gravity  cannot  choose  negative  values  for  gx  or  gy.  Particle  chooses  brx  =  bry  =  1 .  (If 
Gravity  wasn’t  restricted  in  this  way,  it  would  be  quite  a  different  story,  and  Particle  would 
be  advised  to  stick  with  6rx  =  Sry  =  0.)  As  it  is  though,  if  Particle  picks  brx  =  Sry  =  1  then 
no  matter  what  Gravity  chooses,  L  =  1  and  Particle  wins,  getting  to  move  downhill.  (As 
we  said,  this  was  a  pretty  lopsided  game.) 

Suppose  however  that  Particle  is  in  contact  with  an  obstacle  and  must  satisfy  the  motion 
constraint 

brx  +  Sry  <  0  (26) 

along  with  the  regular  bounds  \6rx\  <  1  and  |<5  *>!  <  i-  What  does  Gravity  do  in  this  case? 

Gravity  reasons  as  follows:  “No  matter  what  direction  I  pick.  Particle  always  wants 
to  go  in  that  direction.  Well,  I’ll  fool  Particle  this  time — I’ll  just  pick  a  direction  Particle 
can’t  go  in.  I’ll  choose  gx  =  gy  =  1/2.  That  way,  I  can  guarantee  that  L  <  0 — because  if 
Particle  is  foolish  enough  to  make  6rx  >  0  and  thus  6ry  <  0,  I’ll  switch  to  gy  =  1.  And 
likewise,  if  Particle  makes  6rx  <  0,  I’ll  make  gx  =  1.  But  wait!  What  if  after  I  pick  gx=  1, 
Particle  picks  a  new  direction  so  that  L  >  0?  Well,  I  can  at  least  guarantee  that  Particle 
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Figure  4:  The  motion  of  the  particle  is  constrained  because  it  is  in  contact  with  a  wall.  The 
optimum  strategies  are  to  choose  gx  =  gy  =  1  /2,  and  6rx  =  8ry  =  0. 


can’t  move  downhill  by  sticking  with  gx  =  gy  =  1/2.  There’s  no  way  Particle  can  move  in 
that  direction  at  all.” 

For  its  part.  Particle  reasons:  “If  I  choose  8rx  <  8ryy  then  Gravity  will  pick  gx  =  1 ,  to 
make  L  as  small  as  possible.  But  if  I  choose  8ry  <  8rx,  then  Gravity  will  pick  gy  —  1 .  Aha! 

I’d  better  make  it  so  that  the  smaller  of  the  8rx  and  8ry  is  as  big  as  possible.  But  wait!  If 

I  choose  6rx  >  0,  then  I  have  to  pick  8ry  <  0  because  of  my  motion  constraint.  And  if  I 
choose  8ry  >  0,  then  I’d  have  to  make  8rx  <  0.  Either  way  then.  I’d  be  moving  uphill, 
which  means  L  would  be  negative.  Forget  it!  I’ll  just  choose  8rx  =  8ry  =  0  and  not  move 
at  all!” 

Particle’s  strategy  is  particularly  easy  to  formulate  as  a  linear  program.  Knowing  that 
Gravity  will  focus  attention  on  the  smaller  of  8rx  or  8ry  (and  thus  achieving  a  value  of  L 
equal  to  the  smaller  of  the  two)  Particle’s  strategy  is  simply  to  choose  8rx  and  8ry  such  that 
the  smaller  of  the  two  is  as  large  as  possible  (given  the  constraint  8rx  +  6ry  <  0). 

Now  let  us  apply  this  reasoning  to  equation  (24),  but  with  the  minus  sign  removed: 

max  I  min  grM<5p  I  .  (27) 

j«p>o  1  g>o 

Moo  <  i  \Eg,  =  l  ) 

Consider  the  vector  M<5p,  which  has  length  three  (assuming  we  are  working  with  a  three- 
dimensional  assembly).  Given  a  vector  <5p,  the  g  that  minimizes  grM6p  will  be  such  that 
gi  =  1  where  (M6p),  is  the  smallest  element  of  the  vector  M6p.  Clearly  then,  the  maximum 
of  equation  (27)  occurs  when  6p  is  chosen  so  that  the  minimum  component  of  M<5p  is 
maximized.  Such  a  <5p  can  be  found  by  linear  programming.  Then,  assuming  equation  (24) 
is  true,  the  solution  to  the  dual  of  this  linear  program  will  yield  the  choice  of  g  which 
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maximizes  z(g).  In  the  next  section,  we  explicitly  describe  a  dual  pair  of  linear  programs 
that  enable  us  to  maximize  z( g)  and  prove  that  equation  (24)  holds. 


6  Linear  Programming  Solutions  of  the  Maximin  Problem 

In  this  section,  we  exhibit  dual  linear  programs  to  solve  equations  (21)  and  (23),  but  with 
the  minus  sign  removed.  That  is,  we  will  find  a  solution  vector  g  to  the  problem 


min 

g>0 

Eft-1 

and  a  solution  vector  <5p  to  the  problem 


max  |  mm  l\  izyi 

j«p>o  l  g>o  I 

IM|oo<l  \£ft=i  / 

where  L  =  g7M<5p.  The  goal  is  to  make  the  optimal  solutions  to  the  dual  linear  programs 
satisfy  a  condition  called  a  “saddle-point  condition.”  A  pair  of  vectors  g*  and  <5p*  satisfying 
the  constraints  on  g  and  6p  in  problems  (28)  and  (29)  is  called  a  saddle-point  if  for  all  g  and 
bp  which  also  satisfy  the  constraints  of  the  problems,  the  relation 

g*rM<5p  <  g*rM<5p*  <  grM6p* 


max 

J6p>0 

J*P||oo<1 


(28) 


/  r\r\\ 


holds.  In  terms  of  the  game  of  the  previous  section,  a  saddle-point  indicates  a  pair  of 
strategies  such  that  neither  player  is  inclined  to  change  their  strategy,  providing  the  other 
player  holds  constant  as  well.  We  will  show  that  the  solutions  to  the  dual  linear  programs 
satisfy  the  saddle-point  condition,  and  then  prove  that  any  vectors  that  satisfy  the  saddle- 
point  condition  solve  problems  (28)  and  (29). 

Let  us  define  the  vector  b  by  b  =  (1, 1,  l)7.  In  what  follows,  a  feasible  vector  g  is  a 
vector  satisfying  g  >  0  and  J2gi  —  brg  =  1.  To  express  feasibility  for  a  vector  <5p,  let  I 
denote  the  6 n  x  6 n  identity  matrix,  and  let  e  be  a  vector  of  length  6 n,  with  every  element 
equal  to  one.  If  we  define  the  vector  d  of  length  m+  \2n  and  the  ( m  +  12 n)  x  6/i  matrix  B 
by 

{  0  \  /  -J  \ 


d  = 


and  B  =  |  I 

-I  ) 


(30) 


then  B6p  <  d  implies 


B£p 


-J<5p 

<5p 

-<5p 


< 


(31) 


This  in  turn  implies  J5p  >  0  and  —  e  <  <5p  <  e  which  is  equivalent  to  ||<5p||oo  ^  L  Thus, 
we  will  say  that  a  vector  bp  is  feasible  if  B6p  <  d. 
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6.1  Primal  and  Dual  Linear  Programs 

Let  v  be  a  scalar  and  consider  the  (primal)  linear  program 

maxv  subject  to  (|j  ~  b)  (  tip  )  £  (  d  )  ’  <32> 

If  we  view  the  constraint  bv  —  M<5p  <  0  as  M<5p  >  bv,  we  see  that  to  maximize  v  we  need  to 
make  M<5p  as  large  as  possible.  In  particular,  v  is  bounded  by  the  maximum  that  the  smallest 
element  of  M<5p  can  attain,  given  the  constraints  on  <5p.  Thus,  this  linear  program  exactly 
captures  the  strategy  articulated  at  the  end  of  section  5.2  to  choose  <5p.  Since  setting  v  =  0 
and  6p  =  0  satisfies  the  conditions  in  problem  (32),  it  is  clear  that  a  solution  to  problem  (32) 
exists.  Note  that  for  any  pair  (v,  6p)  which  satisfies  the  conditions  in  problem  (32),  B6p  <  d 
which  implies  that  <5p  is  feasible,  as  defined  above. 

The  dual  linear  program  to  problem  (32)  is 

min  (0,  dT)  ^  ®  j  subject  to  (gT,  sT)  (  J  ~^ )  =  ( 1 , 0),  g,  s  >  0  (33) 

where  s  €  Rm+,2n.  Note  that  any  g  satisfying  the  conditions  in  problem  (33)  satisfies  g  >  0 
and  grb  =  1,  and  is  thus  feasible. 

6.2  Saddle-Point  Condition 

Let  g*  and  s*  be  optimal  solutions  to  problem  (33)  and  let  v*  and  6p*  be  optimal  solutions 
to  problem  (32).  The  duality  theory  of  linear  programming[7]  states  that  the  optimal  values 


of  the  two  linear  programs  are  equal:  that  is, 

V-  =  (0,dr)  (  J."  )  =  dV. 

(34) 

The  duality  theory  also  states  that 

(35) 

The  fact  that  both  vectors  in  this  equation  are  nonnegative  (due  to  the  conditions  in  the 
linear  programs)  means  that 

g*T(M6p*  ~  bv*)  =  0  (36) 

or 

g*rM6p*  =  g*rbv*  =  v* 

since  g*rb  =  1. 

(37) 

Since  g*  and  s*  satisfy  the  conditions  in  problem  (33),  -g*rM  +  s*rB  =  0  or  equiva¬ 
lently,  g*rM  =  s*rB.  For  any  feasible  <5p  then, 


(g*rM)6p  =  (s*rB)5p  =  s*r(B6p)  <  s*rd  (38) 
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since  B£p  <  d  and  s*  is  nonnegative.  Using  the  fact  that  s’rd  =  v*  =  g*rM6p*,  we 
obtain 

g*rM<5p  <  g*rM6p*  (39) 

for  all  feasible  6p. 

To  obtain  the  other  half  of  the  saddle-point  condition,  observe  that  M6p*  >  bv*.  Since 
g  >  0  and  grb  =  1  for  any  feasible  g, 

gr(M£p*)  >  gr(bv*)  =  (grb)v*  =  v*  =  g*rM<5p*.  (40) 

Combining  the  previous  two  inequalities  gives  the  saddle-point  condition 

g*rM<5p  <  g*rM(5p’  <  grM<5p*.  (41) 


for  any  feasible  g  and  6p. 

6.3  Maximin  Result 

We  claim  that  problem  (28)  attains  its  minimum  value  (and  thus  z(g)  attains  its  maximum 
value)  for  g  =  g*,  and  that  problem  (29)  attains  its  maximum  value  for  6p  =  6p*. 
Furthermore,  these  values  are  equal.  In  what  follows,  we  restrict  our  attention  to  feasible 
vectors  6p  and  g  (that  is,  we  will  not  write  out  the  feasibility  conditions). 

To  show  that  g*  minimizes  problem  (28)  we  need  to  show  that 

max  g*rM<5p  <  max  grM<5p.  (42) 

ip  ip 

for  all  feasible  g.  From  the  saddle-point  condition,  g*rM6p  <  g,7M6p*  for  all  feasible  <5p, 

so  max  g‘rM<5p  is  bounded  from  above  by  g*rM6p*.  Since  g*rM<5p*  <  grM6p*  for  any 
ip 

feasible  g,  we  have 

max  g*rM(5p  <  g*rM<5p*  <  grM6p*  <  max  grM6p 

ip  ip 

for  any  feasible  g.2  Thus, 

min  ( max  grM6p )  =  max  g*  rM<5p  <  g*rM6p*. 

g  \  ip  /  ip 

Similarly,  the  saddle-point  condition  yields 

nun  grM<5p  <  g*rM6p  <  g*rM(5p*  <  mingrM6p* 

for  all  feasible  <5p,  so  <5p*  maximizes  problem  (29)  and 

g*rM6p*  <  mingrM£p*  =  max  (mingrM6p  j . 

g  ip  \  *  ) 

2The  inequality  grM6p*  <  maxip  gTM6p  follows  from  the  fact  that  max$p  grM6p  is  greater  than  or  equal 
to  grM5p  for  any  choice  of  <5p.  In  particular  then  if  we  choose  <5p  =  <5p*,  maxip  grM£p  >  grMrfp*. 


(43) 

(44) 

(45) 

(46) 
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This  proves  that  g  =  g*  and  6p  =  <5p*  are  optimal  for  problems  (28)  and  (29)  respectively, 
and 

min  (maxgrM<5p )  <  g*rM<5p*  <  max  (  mingrM6p)  .  (47) 

*  \  «P  /  «P  \  8  J 

Thus,  we  can  find  a  value  of  g  which  maximizes  z(g)  by  solving  the  linear  program  of 
equation  (33).  The  motion  6p  =  <5p*  which  minimizes  6U  for  the  orientation  g*  according 
to  the  constraints  ||^p||oo  <  1  and  J6p  >  0  can  be  found  by  solving  the  linear  program  (32). 

Finally,  we  can  use  the  optimal  vectors  g*  and  <5p*  to  obtain 


max  (  min  grM<5p )  =  min  grM<5p*  <  g*rM<5p*  <  max  g*rM<5p  =  min  ( max  grM6p ) 

£p  \  g  J  8  6p  g  \  6p  / 

which  yields 

max  (  minL  1  <  min  (  maxL )  .  (48) 

«p  V  g  J  ~  g  \  Sp  J 

Combining  the  inequalities  in  (47)  and(48)  proves  that 


min 

g 


^maxgrM6pj  =  g*rM6p*  =  max  ^mjn  g 


(49) 
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